Non classicality of generalized vortex states produced by multi-photon subtraction 

from two mode squeezed vacuum 
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We discuss the non classicality of a class of vortex states which are engineered by photon substrac- 
tion from two mode squeezed vacuum. We calculate the log negativity parameter, the von Neumann 
entropy of the reduced density matrix and the concurrence. The negativity of the Wigner function 
is discussed as a measure of the non classicality and the quantum interference pattern obtained 
therein as a possible measure of the entanglement present between the two modes. This measure of 
entanglement is compared with the results obtained from concurrence. 
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I. INTRODUCTION 



Over the last decade, the physical characteristics of the 
entanglement of quantum mechanical states, both pure 
and mixed, has been recognised as a central resource in 
various aspects of quantum information processing pQ. 
Given an arbitrary bipartite pure state \^ab), the en- 
tropy of entanglement E{\^> ab)) 0j namely the von 
Neumann entropy of the reduced density matrix pa = 
Tr b\^ab){^ab\, tells us exhaustively about the prob- 
abilities of transforming I^ab) using LOCC, into other 
pure states in an asymptotic sense. Many efforts have 
also been devoted to the study of mixed state entangle- 
ment [SI H]. Given the central status of entanglement, 
the task of quantifying the degree to which a state is en- 
tangled is important for quantum information processing 
and correspondingly several measures such as entangle- 
ment of formation [3, 5], von Neumann entropy [6] and 
negativity [TJ [7] have been formulated. The entangle- 
ment of formation or concurrence have been the topic of 
much interest lately. Although originally formulated to 
quantify pure state entanglement, it was later extended 
to mixed states [8HTT]. 

Among all (quasi) probability distributions in quantum 
mechanics [H], Wigner function [13 stands out because 
it is real, non singular, yields correct quantum mechanical 
operator averages in terms of phase space integrals and 
possesses definite marginal distribution [T?]. Moreover, 
the negativity of Wigner function depicts non classicality 
of the state [15] , which we highlight in this article as a 
measure of entanglement also. 

Optical vortices or phase singularities in the field of light, 
described by a non-separable two-dimensional field, have 
drawn a great deal of attention in the last two decades 
[16) and these have prompted to start a new branch in 
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physical optics known as singular optics [T7]- Because of 
their specific spatial structure and associated orbital an- 
gular momentum (0AM), they find various applications 
in the field of optical manipulation |18| , optical communi- 
cation |19j . quantum information and computation |20| . 
For a vortex, the OAM follows the same numbers as vor- 
ticity or topological charge of the vortex i.e. an optical 
vortex of topological charge m, carries an OAM of mh 
per photon. Most applications of vortices lie in the clas- 
sical realm. However, in recent literature [21H25J quan- 
tum optical vortex has also been studied. It has been 
shown that such states can be generated from two mode 
squeezed vacuum under a linear transformation belong- 
ing to the SU(2) group with some restrictions. 
We studied the Wigner function of a general elliptical vor- 
tex, the circular vortex being a special case of an elliptical 
vortex which is constructed from two separate squeezed 
vacuum states by mixing them through a beam split- 
ter (BS) or a dual channel directional coupler (DCDC) 
[2"3"] . We also studied the general elliptical vortex from 
an information theoretic viewpoint [2B]. We looked at 
the suitability of it for carrying information by studying 
the von Neumann entropy. 

Recently Agarwal [23] analyzed quantum optical vortex 
by subtracting a photon by a 99% transmitting beam - 
splitter (BS) from one of the two modes of spontaneous 
parametric down-conversion (SPDC), generated by the 
non - linear crystal (NLC). He showed that the subtrac- 
tion from one mode results into an addition to the other 
mode; it produces a vortex, which is elliptical in shape. 
It must be noted that a single photon subtraction or ad- 
dition has already been demonstrated by Parigi et al. 
[27] . A more recent paper studied the enhancement of 
entanglement, though defined differently, in photon sub- 
traction/addition [55] 

The paper is organized as follows. In section [H] we in- 
troduce the generalized vortex state arising from multi 
- photon subtraction of the two mode squeezed vacuum. 
We show that multi - photon subtraction indeed gives rise 
to vortex states of higher order. We discuss the non clas- 
sicality of this state with the help of the Wigner function. 



In section III we discuss the entanglement properties of 
the generalized vortex state with the logarithmic nega- 
tivity and entanglement of formation. We discuss the 
variation in mixedness with squeezing. We also present 
an interpretation of the quantum interference patterns 
obtained from the Wigner function as a measure of the 
entanglement present between the two modes by com- 
paring it with the results obtained from our study of 
concurrence. We conclude this article after pointing out 
the significant results and directions for future work in 
section IIVI 



mode a (signal) and y, p y correspond to mode b (idler). 
The quadrature distribution thus obtained has the fol- 
lowing form 
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The non-separability is evident from the fact that Eq 



II. VORTEX STATES ARISING FROM 
PHOTON SUBTRACTION/ADDITION 

Quantum optical vortices are phase space singularities 
in the quantized radiation field which possess orbital an- 
gular momentum. There are methods of producing these 
states. One of them is by mixing two squeezed vacuum 
modes through BS or dual channel directional coupler 
[2"5] . The other one is by photon subtraction using a BS 
from one of the outputs of a SPDC [24l [27] as shown in 
Fig. [T] Vortex of higher orders can be produced by in- 
creasing the number of BSs and avalanche photo diodes. 
Here we discuss the properties of vortices created by sub- 
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FIG. 1: A schematic of the experimental setup for the pro- 
duction of vortex state of order 2 by subtracting two photons 
from a two - mode squeezed vacuum using two 99% beam 
splitters(BS 1 and BS 2). NLC is a nonlinear crystal. SPDC 
generates the signal (mode a) and the idler (mode b). APD 
1 and 2 are Avalanche Photo Diodes. 

tracting photons. 

A two-mode squeezed beam is represented by 



10 = exp (fcW - £* a6)|0,0), £ = re 1 * 



(1) 



where £ is a complex parameter, r is the squeezing pa- 
rameter and a and b are the regular bosonic mode op- 
erators satisfying the commutation relations [a, a'] — 
[b, tf] = 1, [a, tf] = 0, etc. The corresponding quadrature 
distribution is obtained by replacing the bosonic mode 
operators by the quadrature operators x,y,p x and p y 
which obey the relations [a;,^] = [y,p y ] = i- Here we 
consider H = 1. The quadratures x, p x correspond to 




(c) 

FIG. 2: (Color online) Plot of the intensity (a) 3D plot, (b) 
contour plot and (c) phase oi^%(x,y) for fc=3. The singular- 
ities are clearly seen in Fig. (c). 



(pi) cannot be written as a product of two functions each 
of which are solely dependent on either x or y. 
Now one photon can be subtracted from mode b with 
the help of a beam splitter. 



\0 is) = Abexp (fcW - Cab) |0, 0) 



(3) 



where A is the constant of normalization. This photon 
can be detected using a single photon detector. The prob- 
ability of detection is proportional to the square of reflec- 
tivity |24j . If an arbitrary number of photons, say k, are 
subtracted from the two mode squeezed vacuum, Eq. ffib 
can be generalized to the following 

\0{ S) = A k b k exp (frt&t - Cab) |0, 0) (4) 

We can rewrite Eq. Q as follows 

= A k b k - 1 exp i&tf - Cab) 
x exp (-^aW + Cab) b exp (faW - Cab) |0, 0) 

(5) 



\ci s) 



This can be simplified by repeated use of the two mode 
Bogoliubov transformation to finally arrive at the follow- 
ing simplified form 
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A k exp (^a t 6 t - Cab) (bcoshr + a t sinhre i9 ) fc |0, 0) 
Eq. ([6]) can be further simplified to 

\C j>> = A k exp (frW - Cab) \k, 0) sinh fc re lke (7) 



One can also write the state (Eq. (|7j)) in terms of Fock 
states by using the decomposition of the squeezing oper- 
ator and obtain 



i \ e lke 
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cosn r 



(8) 



So we see that the state ((8) can also be obtained by 
adding k photons to |£) in the mode a. This is similar 
to Parigi et aVs experiment in which they added a single 
photon to generate the vortex state. The state ^ ex- 
hibits the property that there is a fixed difference in the 
number of photons between the signal and idler states. 
In literature this is one of the so called "pair coherent 
states" [351 [3U] • Pah coherent states provide an impor- 
tant example of non classical states of the two mode ra- 
diation field. It has been studied in detail for their non 
classical properties and as examples of EPR states [29] . 
The other property of interest is the vortex structure. 
The quadrature distribution for the state ^ can be ob- 
tained by writing the bosonic creation operator a^ in 
terms of x a and d/dx a as follows 



* k s) (x,y) 



ihii 



2 fc / 2 coshr 



d_ 
dx 



¥ s \x,y) (9) 
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FIG. 3: (Color online) Plot of the intensity (a) 3D plot, (b) 
contour plot and (c) phase of ^(as, y) for k=4. The singular- 
ities are clearly seen in Fig. (c). 



where ^^ s '(x, y) is the quadrature distribution of the two 
mode squeezed vacuum, Eq. ([2]). For k = 1, Eq. Q can 
be expressed as follows 

*i {x,y) = 



(1 — if) ' 7r 1 / 2 cosh r 



xexp 



2xyy - (x 2 + y 2 ) rj 2 _ 1 2 

1 - rf 2 [ +V ' 



(10) 



This is seen to possess a vortex structure for rj = i\r]\. 
In other words, the method of single photon subtraction 
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FIG. 4: Wigner function of the TMSV state, for the five photon subtracted case (k=5). 



from one of the two modes taken from the output of a 
spontaneous parametric downconverted beam results in 
the formation of a vortex state of order 1. Intuitively 
it can be argued that the process of subtracting k pho- 
tons should give rise to vortex states of higher order. As 
has already been pointed out, multi-photon subtraction 
from one of the modes of the two mode squeezed vacuum 
is similar to adding an equal number of photons in the 
other mode. The basic idea is to have a difference in the 
photon number between the two modes. This gives rise 
to higher order vortex states where the order or vorticity 
is determined by the difference in number of photons in 
the two modes. We call these states Two Mode Squeezed 
Vortex (TMSV) states. 

We study the states arising from 3 photon and 4 photon 
subtraction in Fig. ^ and Fig. ^ respectively. It is 
evident that this method gives rise to higher order vortex 
states. The order can be determined from the phase plots 
by counting the number of singular points. For k—l, the 
results of 24J is reproduced. It is observed, with increas- 
ing order, the intensity distribution associated with the 
vortex begins to divide into separate lobes, as can be seen 
from Fig. (pL when three photons are subtracted (fc=3). 
When another photon is subtracted, the corresponding 
intensity distribution is completely divided into two dis- 
tinct lobes. These results serve as a pointer to study the 
entanglement properties of these states and its variation 
with the order of the vortex. 

The non classical nature of these states can be studied 
with the help of the Wigner function. For a number state 
|n), it is written as |31j . 

W(x,y)= 2 (-l) n £ n (U 2 )e- 2 > 2 , (11) 

where 4? = x 2 + y 2 and C n is the Laguerre polynomial 
of order n. The Wigner function for a two mode Fock 



state, \n,m), can then be obtained in terms of coherent 
state parameters a — x — ip x and /3 — y — ip v as follows 



W{a,P) 



-l) n £„[4|a| 2 ]exp(-2; 



a 



x (-l) m £ m [4|/3| 2 ]exp(-2|/3| 2 ) (12) 



The Wigner function for the two mode state |n,0) can 
then be simplified and written as 

W(a, 0) = -1 (-1)" C n [4|a| 2 ] exp [-2 (|a| 2 + |^| 2 )] . 

It has been shown that if two density matrices are related 
by a squeezing transformation, then the Wigner functions 
are related by the following transformation [24] 



W(a,P) = W(aJ), 



where, 



coshr 
— sinhr e 1 ' 



-sinhr e 10 
coshr 



(14) 



(15) 



So, the Wigner function for the TMSV, |£)jj. , can be 



written as, 

W(a, 0) = ~ (-l) k C k [4|5| 2 ] exp [-2 (|<5| 2 + |£| 

(16) 
where k is the number of photons subtracted. In Fig. Q 
we study the variation of the Wigner function as a func- 
tion of two of the quadratures out of four (two quadra- 
tures for each of the two modes), keeping the other two 
quadratures as zero. 

The non classicality of the TMSV states is evident from 
the nature of the Wigner function which shows negative 
regions for all the photon subtracted states. In Fig. (4a) 



and Fig. (4b), we study the cross correlation between 



the same quadratures of two different modes. It is ob- 
served, though not shown in the figure, that the peak 
alternates between 0.4 (for even k) and -0.4 (for odd k). 
It attains the positive maximum when an even number of 
photons are subtracted. Negative maximum is attained 
when an odd number of photons are subtracted. The vol- 
ume of the negative region of the Wigner functions can 
be further studied as a measure of the non classicality 
of these states [15] . It can also be interpreted as a mea- 
sure of entanglement between the two modes as explained 
in section IIIC Fig. (4c) and Fig. (4dl represents the 
cross correlation between two quadratures of the same 
mode. These are exactly similar in nature to that of the 
cross correlation between same quadratures of different 
modes. In Fig. (4e) and Fig. (BfJ, we study the cross 



correlation between different quadratures of two different 
modes. These show interesting quantum interference ef- 
fects. Similar interference pattern has been reported in 
[2"2l 152] . We present an interpretation of this quantum 
inteference later in this article. 
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(a) Two photon subtraction 
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III. ENTANGLEMENT PROPERTIES 

In this section we study the mixedness of the TMSV 
states as well as the entanglement properties of them. 

A. Logarithmic negativity 

There exists various measures of entanglement. A 
much used tool is the log negativity parameter pQ. We 
use this to quantatively study the entanglement in the 
vortex state. It is defined by 



log 2 (l + 2A0 



(17) 



where Af is the modulus of the sum of all the negative 
eigenvalues associated with the partial transpose p PT of 
the density matrix |£)j. ^ (CI- The density matrix of the 
TMSV state can be written from Eq. ^ as follows 

P = H c mc n \m + k, m)(n + k, n\e li - m ' n)e (18) 



where, 



tanh p : 



cosh r 



(19) 



Here p = m, n. Therefore, the log negativity parameter 
of Eq. (17) becomes 



£fc 



l°g 2 (H c ») 



(20) 



(b) Three photon subtraction 
FIG. 5: Entanglement ratio tk for the TMSV states 



uum gives the result log 2 (e 2r ) using c n = tanh"r/coshr. 
The entanglement in the vortex state can be compared 
with that in the state |£) by studying the ratio of log 
negativities of the two states. The ratio has the form 



Cfc 



Z^' 



(21) 



We study the results in Fig. ^. For k = 1, i.e. single 
photon subtraction, the results of [23] were reproduced. 
It was observed that for two photon subtraction there 
was a very small rise in entanglement which falls sharply 
after reaching the peak. Further subtraction worsens 
the situation. For three photon subtraction the initial 
value falls to half. For four photon subtraction, the 
starting value is 0.042. This seems pretty interesting 
given the nature of the higher order vortex states in Fig. 
^ and Fig. ([3|. When three photons are subtracted, 
it gives rise to a vortex state of order 3 as can be seen 
from Fig. pc|. But it is also evident from Fig. (2b) 



that the corresponding intensity distribution almost 
disintegrates into separate lobes. The situation worsens 
when four photons are subtracted (Fig. (3c) and Fig. 
(3b)). The intensity distribution, in this case, is seen 



to have completely divided into two distinctly separate 
lobes. Although the lobes are separate, it does not imply 
that there is no entanglement between photons sampled 
from the different lobes. 



A similar calculation for the two mode squeezed vac- 



B. von Neumann entropy 



C. Entanglement of formation 



The von Neumann entropy, also known as the entropy 
of entanglement is defined as 

S a = -Tr (p a log 2 Pa) = - ^2 ^n log 2 A„, (22) 



where, A„ are the eigenvalues of p a , the reduced density 
matrix of |£)j. . The reduced density matrix is defined as 

Tr a [|Oi S) i S) (CI 
as follows 



Using Eq. (18), we can write it down 
p„ = 53cg|m + l)<m + l| > (23) 



where Cfc,s are the same coefficients as in Eq. (19). The 



entropy of entanglement, as is well known, can be in- 
terpreted as a measure of the mixedness of the states. 
Its values, when appropriately scaled, vary between zero, 
for pure states, and one, for completely mixed states. We 
study the variation of the same with the squeezing pa- 
rameter, in Fig. ^. 
We observe, that the entropy increases with the squeez- 
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FIG. 6: Entropy of entanglement for the TMSV states for 
fe=4. 



ing parameter, r. It increases from to reach a max- 
imum value, at r = 3.2. Since von Neumann entropy 
is interpreted as a measure of the mixedness [10] . this 
implies that the states change from completely pure to 
completely mixed state with the increase in squeezing. 
When r — 0, that is, there is no squeezing, the states 
remain pure. It gradually looses its purity with increas- 
ing r and becomes completely mixed at a higher value 
of r. The dependence of the mixedness of the TMSV on 
the squeezing parameter is not very surprising. In our 
setup, we use the SPDC to squeeze the two mode vac- 
uum. This also results in the formation of a mixed state. 
So, the squeezing parameter apart from parametrizing 
the squeezing of the two mode vacuum also parametrizes 
the mixing of the two modes. It is also observed that the 
variation of the von Neumann entropy with squeezing 
parameter is independent of the order of the vortex. 



In this section we study the entanglement of the TMSV 
state with the help of concurrence 012], also expressed as 
the entanglement of formation [3] ■ It is a measure of the 
amount of entanglement needed to create the entangled 
state. It is defined as jTD] 

C (p) = min Pi ^ Y,Pi E (l^:)(^l) (24) 

i 

The minimization is taken over all possible decomposi- 




FIG. 7: (Color Online) Entanglement of formation of the pho- 
ton added (subtracted) vortex state, m denotes the number 
of photons added (subtracted) in one of the modes. 

tions of pure states ipi with probabilities pi, which taken 
together reproduce the density matrix p of the given 
mixed state. On the other hand, E ({ipi) (ipi\) is the en- 
tropy of entanglement, more commonly known as the von 
Neumann entropy, of the reduced density matrix as de- 
fined in the previous section IIIB Using Eq. (18), Eq. 



( 23 ) and Eq. ( 22 ) , we arrive at the following simplified 



form of the entanglement of formation 



C (P)k = min r 



• IOi 



n,m 



X»S: 



2 c p 



(25) 



where, q's are the same as in Eq. (19). We study 



the variation of the entanglement of formation with the 
squeezing parameter and the order of the vortex in Fig. 
ffh. We observe that the C(p) k increases with increase 
in r, the squeezing parameter. The rate of increase grad- 
ually slows down till it reaches the maximum value and 
then falls of rapidly to 0. It should be noticed that C {p) k 
attains its maximum value at r ~ 2.1. This means that 
the two modes become more and more entangled with 
increasing r until it reaches a maximum value beyond 
which the modes start becoming disentangled. We in- 
fer that there exists an optimum value of the squeezing 
parameter which when applied at the beam splitter pro- 
duces maximum entanglement of the two modes. If we 
go on squeezing the modes more and more, they will no 
more be entangled. As we move to higher orders of the 
vortex, the range of values of r which produces maximal 
entanglement, starts broadening. At k=5, there exists 



W(x,Py?- 
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(c) r=2.1 
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(d) r=3 

FIG. 8: Wigner function W(x,p y ) of the TMSV states show- 
ing quantum interference for different values of r, the squeez- 
ing parameter, with k — 5. 



a well defined range of r that produces maximal entan- 
glement of the two modes. The entanglement falls off 
rapidly on both the ends of this range. 
We present here an interpretation of the quantum in- 
terference effects observed in section |n} In Fig. ([8]), we 
study the effect of the squeezing parameter on quantum 
interference. It is observed that the interference is most 
prominent for r = 2.1 (see Fig. 



8c). This observation 



has a striking similarity with the results obtained from 
our study of the entanglement of formation. The inter- 
ference effects decrease as we move away from r = 2.1 on 
either side which is exactly similar to what we obtained 
from the entanglement of formation. Broadening of the 
range of r which gives rise to a well defined interference 
fringe system with increase in the order of the vortex was 
also observed. Hence, we suggest that the quantum in- 
terefence effects arising from the Wigner function of the 
cross correlation between different quadratures of the two 



modes can be interpreted as a signature of the entangle- 
ment present between the two modes. It is seen that the 
number of fringes as well as the volume of the fringes 
change with changing r. However, the maximum num- 
ber of negative fringes is restricted by the order of the 
vortex. Since the negativity of the Wigner function is 
interpreted as a measure of non classicality, we suggest, 
the volume of the negative part of the quantum interfer- 
ence fringes might be a possible candidate for quantifying 
entanglement and needs to be further studied. 



IV. CONCLUSION 

In conclusion, we have studied the non classical proper- 
ties of the photon subtracted (added) two mode squeezed 
vacuum state. We have shown that the order of the vor- 
tex created by photon subtraction (addition) is equal to 
the number of photons subtracted (added) from one of 
the modes. Rather it should be stated that the order was 
equal to the difference in the number of photons between 
the two modes achieved by subtracting (adding) photons 
from one of the two modes. The intensity distribution of 
the higher order vortex, thus created, starts disintegrat- 
ing. This gives rise to two distinctly separate lobes. 
The non classicality was evident from the negative re- 
gions of the Wigner function. The volume of the neg- 
ative region can be interpreted as a measure of the non 
classicality and should be a topic of further investigation. 
We studied the mixedness of these states with the help of 
von Neumann entropy and the entanglement properties 
using the logarithmic negativity and the entanglement 
of formation. The states were observed to change from 
pure to mixed states with increase in squeezing. It was 
found that there exists a threshold value of the squeezing 
parameter, beyond which the states become completely 
mixed. It was observed that the entanglement decreases 
with increasing order of the vortex with the help of the 
logarithmic negativity. However, no such evidence was 
found while studying the entanglement of formation or 
concurrence. On the contrary, it was observed that with 
increasing order of the vortex, the range of values of r, the 
squeezing parameter, which gives rise to optimum entan- 
glement, broadened. It was found from the study of the 
entanglement of formation that the maximum entangle- 
ment achieved for all orders of the vortex was equal. The 
only difference was the range of values of r for which it 
was possible to achieve the maximum entanglement. The 
higher the order of the vortex, the broader was the range 
of r. 

It is also suggested that the quantum interference effects 
obtained from the study of the cross correlation between 
different quadratures of two different modes might be in- 
terpreted as arising due to entanglement between the two 
modes. The volume of the negative regions of the inter- 
ference fringes need to be further investigated as a pos- 
sible measure of the entanglement. This further points 
to the need of a more robust measure of entanglement. 



The present measures of entanglement are found to be 
inadequate in that they are dependent on the system of 
states being studied and are not universal measures. The 
TMSV state discussed in this article should also be fur- 
ther investigated to see how it behaves under teleporta- 
tion. 
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